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Abstract 

We develop a notion of computability and complexity of functions over the reals, 
which seems to be very natural when one tries to determine just how "difficult" a certain 
function is. This notion can be viewed as an extension of both BSS computability 
|B(ySS98] and bit-computability in the tradition of computable analysis IWeiOOj as it 
relies on the latter but allows some discontinuities and multiple values. 

1 Introduction 

The main goal of this paper is to provide a simple definition which would capture the intuitive 
notion of whether / : M ^ R is an "easily" computable function. 

There are two main currently existing approaches to the computability and complexity 
of real functions. One is the BSS approach, where algebraic operations are performed on 
real numbers that are stored with infinite precision. This approach is described in 



and |B(ySS98j . The other approach, which we call bit-computability, goes all the way back 
to the Polish school in the 1930-50s, and can be formulated best as: "given a good rational 
approximation of x, compute a good rational approximation of /(x)" . More recent references 
on the subject include |Ko91j and |WeiOOj . 

The quality of a given definition for computability of real functions (or any notion of 
computability in general) can be judged by how well it matches the intuitive notion of 
"easy" vs. "hard" vs. "impossible" . For example, in the discrete case, it it much easier to 
add two numbers x + y (an operation that takes linear time in the size of the input), than 
factor an integer n (an operation most believe requires time superpolynomial in the size of 
the input), and solving the Halting Problem is truly impossible (for example, the Goldbach 
conjecture can be presented as a simple instance of the Halting Problem). We keep this in 
mind while examining the different approaches to the computability of real functions. 

We consider some reservations to both approaches mentioned above. One complaint 
about the BSS approach is that it is too focused on the algebraic simplicity of the function 
/, allowing only piecewise-semialgebraic functions to be computable. In particular, simple 

"'^ Research is partially supported by an NSERC postgraduate scholarship 
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functions such as x ^— and x i— > ^/x are not computable in this model (cf. |Brt03bj ). 
There are also problems in the opposite direction: this model classifies all the constant 
functions x ^ a, which are extremely simple algebraically, as computable. However, there 
are infinitely many real a's such that computing a with an arbitrarily good precision would 
allow us to solve the Halting Problem. We will discuss more of these issues, and ways to 'fix' 
them in section El below. 

On the other hand, in the bit computability model all the computable functions must be 
continuous. Thus even the simplest step function so(a;) = 1 if a; > and /(x) = if a; < 
is not computable in this model. This function is extremely simple from the mathematical 
point of view (involving only the constants and 1 and the comparison operations in its 
description), hence it would be reasonable to expect it to be computable at least in some 
sense. For a real function / it is hard to say whether it is uncomputable in this model 
because it is "difficult" in some profound way, or because it is just discontinuous in one or a 
few points. It would also be nice to have a setting in which, by analogy with the discrete case, 
we have a strong connection between the computability of a set A and the computability of 
its characteristic function xa- Another problem is how to properly define the computability 
of multi-valued functions such as ^ : C — C. It is obvious that such a simple function 
should be computable, but it doesn't have a continuous (and hence bit-computable) branch 
defined on the entire complex plane C 

We deliberately restrict our discussion to the simple case where the functions are from 
from a simple rectangle such as [0, l]'^ to M^. The goal being not to give the broadest 
definitions and prove the most general theorems, but rather provide a simpler discussion ac- 
cessible to a broader audience. We will try to use only basic background in the computability 
and complexity theory and in the topology of M*"'. 

The paper is organized as follows. In section |21 we give an introduction to the bit model of 
computation for sets, along with some properties and examples. This section can be viewed as 
a separate simple introduction to the subject. In section El we discuss the BSS computability 
for sets, and propose three natural modifications to the model that make BSS-computability 
equivalent to bit-computability. In sectionlHwe introduce another notion of set-computability 
- weak computability, and show its equivalence to the standard bit-computability. Sections El 
and m can be read independently. In section El we use results from sections El and El to propose 
a new computability and complexity definition for real functions extending the computability 
in both models, and the classical complexity in the case of continuous functions. 

Acknowledgement. The author wishes to thank his graduate supervisor, prof. Stephen 
Cook, for his insights and encouregement during the preparation of this paper, and for the 
countless helpful discussions. 
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2 The Bit Model 



2.1 The Model of Computation 

The comput ability of functions in the bit model as we know it today was first proposed by 
Grzegorczyk |Grz55j and Lacombe |Lac55j . It has been since developed and generalized. 
More recent references on the subject include |Ko91j . |PR89j . and |WeiOOj . 



The basic model of computation here is a Turing Machine. One can think about it as 
a program in any programming language on the computer. We will usually denote Turing 
Machines by the letter M. Obviously, a computer has a finite memory and cannot store a 
whole real number. Instead, a naming system is used to represent a real number x. The 
most popular naming system for R uses the dyadics ro = {^|mGZ, nG N}. A name for a 
number x would be a sequence of dyadics 0(1), 0(2), 0(3), . . . such that |0(n) — x\ < 2^". In 
particular, (j){n) — > x as n ^ oo. Note that the same real number x can have more than one 
name, e.g. 0, |, |, ^, || . . . and I? ^! H • • • are both names for the same number a; = |. 

The oracle terminology is just a natural way to separate the complexity of computing x 
from the complexity of computing on a; as a parameter. For most purposes on can think of 
the oracle for x as an infinite tape containing the binary expansion of x. 

Consider a function / : R — R. In plain language a program M computing / would 
output /(x), provided the input x. More precisely, given a name for x, M should output 
a name ip for f{x). Of course, both the name of x and of f{x) take infinitely long to write 
down in general. We deal with this problem as follows: M is allowed to query (pim), which is 
a 2~™'-approximation of x for any natural m, and is required to output ip{n), where n is given 
to M as a parameter. In other words, M is allowed to get arbitrarily good approximations 
of X to compute f{x) with a given precision 2~". Note that M does not care about the way 
0(m) is obtained or computed, thus we say that is an oracle representing x to M, and we 
write M'^ to emphasize this fact. We sometimes write M'f'{n) to emphasize the fact that M*^ 
gets one precision parameter n. 

The definition extends naturally to a function f : R'' ^ R^ Here M = M'l'^''^^'-^^'={n) is 
allowed to query each of the k parameters with an arbitrarily good precision and is required 
to output the i values of / with precision 2~". 

Example: Compute the function / : [0, 1] — > [0, 1], f{x) = x"^ . 

Solution: On an input n, query for q = (f){n + 1). Output the dyadic number g^. 

To show that it works, we need to see that |g^ — < 2^". We know that jg — x| = 
|0(n + 1) — x| < 2~"'~^, hence 

- = \q + x\ ■ \q - x\ < 2 ■ \q - x\ < 2 ■ 2-"-^ = 2^". 

The running time of the machine M'^{n) is the largest running time over all legitimate 
oracles 0. 
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2.2 Basic Properties and Examples 



One of the main properties of computable functions is that they are continuous. To see 
this assume that / is computable and let x be a point in the domain of /. Let (f) be some 
oracle representing x as described above. The computation M'^{n) terminates after finitely 
many steps with an output q such that \q — f{x)\ < 2~". M only queries (j) with some finite 
precision 2""^, and the computation would be the same for any oracle which agrees with 
(j) on the first m values. If we choose and carefully, we see that \q — f{y)\ < 2~" for all 
y such that \x — y\ < 2^'"^^. This shows that / must be continuous. 

Moreover, one can show that if the domain of / is a simple compact (closed and bounded) 
set, such as the interval [0, 1], then the modulus of continuity of / is computable. That is, 
we can compute a function m(n) such that |/(a;) — f{y)\ < 2""" whenever \x — y\ < 2~"^^^\ 

We recall the definition of computable real numbers, introduced by Turing in |Tur86j . In- 
formally, this definition says that a real number x is computable if we can compute arbitrarily 
good approximations of x. 

Definition 1 A number x E W is computable if and only if a representation (p of x as 
described above can be computed. 

Most "standard" continuous functions are computable in this model. For instance, the 
exponential function f{x) = is computable on R using the Taylor series expansion of 

= ff- hard to estimate the number of terms and the precision of x we need 

to consider in order to get a 2~"-approximation of e^. 

Similarly, polynomial, rational and trigonometric functions with computable coefficients 
are computable (on properly chosen domains). Moreover, if /, g are computable on the same 
domain D and c is a computable constant, then c ■ /, f + g and / ■ g are also computable. 
^ is computable if g ^ on A. 

Let a G -D be a computable number, and /(a) = g{a), then the piecewise defined function 



is also computable. Note that the condition /(a) = g{a) is essential here, for otherwise the 
function h would not be continuous (and cannot be computable in this case). 

It must be noted that in the continuous case, only truly "pathological" functions are 
uncomputable. This notion of computability seems to be effective in classifying continuous 
functions. 

One of the possible disadvantages of this model is the fact that even the simpliest dis- 
continuous function, the step function 



is not computable under this definition. This is in contrary to the intuition that Sq must be 
a "simple" function. One could also argue that some physical systems, e.g. quantum energy 
levels, are best described using step functions and other simple discontinuous functions. One 
of the goals of the present work is to develop notions which deal with this problem. 




X G D, X < a 
X E D, X > a 




X <0 
x> 
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2.3 Complexity of Real Numbers and Real Functions 



The time complexity Ta{n) of a number a is the time complexity of computing a 2~"- 
approximation of a. Note that this is the time the fastest machine would take to compute 
a. If Ta{n) < p{n) for some polynomial p{n), then we say that a is poly-time computable. 

Example: Let's consider the time complexity T^(ra) of computing the number vr. One 
way to compute tt is to use the Taylor series of arctan: 

vr = 4 ■ arctan(l) = 4 ■ ("l - 1 + 1 - 1 + . . = 4 • ^"^^^ 



3 5 7 / ^2i + l 



It is not hard to see that one would need to take exponentially many terms of the series to 
obtain a 2~"'-approximation of vr. On the other hand, one can use the Bailey, Borwein and 
Plouffe jBBP97j formula for vr: 



oo ^ 



1 1 



16* V8z + 1 Si + 4 82 + 5 82 + 6 

i=l 

(or any one of the many other exponentially convergent vr expansions). One needs only 
linearly {k. j) terms of the BBP series to get a 2~"-approximation of vr. Hence vr is, in fact, 
poly-time computable. 

The definition of the complexity of a function / arises naturally from the definition of 
computability above. We define the worst case time it takes to compute f{x) with precision 
2~" to be the time complexity of f{x), and denote it by Tf{n). Note that Tf may depend on 
the domain D of /. In general, the time complexity may increase as the domain expands. 

If Tf{n) < p{n) for some polynomial p{n), we say that the function f{x) is poly-time 
computable on D. The standard functions mentioned in the previous section are typically 
poly-time computable, using the standard numerical-analytic techniques. 

Note that the time complexity of the constant function f{x) = a is equal to the time 
complexity of the number a (we just ignore the argument x). 

Example: The function f{x) = is poly-time computable on the domain [0, 1]. To see 
this, we use the Taylor series expansion 

oo ,• 
X 



i=0 



It is not hard to see that on the [0, 1] interval, n (and even fewer) terms of the series suffice 
to get a 2~'^-approximation of f{x). All we need to get the sum of the first n terms is 
0{n) additions and multiplications on 0(n)-bit numbers. Hence even with the naive 0{n'^) 
multiplication algorithm, we obtain Tf{n) = 0{n^). 



2.4 Computability and Complexity of Real Sets 

According to |BW99j . probably, the first definitions of effective subsets of based on 
the concept of computability have been proposed by Kreisel and Lacombe in 1957 |KL57j . 
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jLac58j . We refer the reader to [BW99j and |WeiOOj for a more detailed discussion. By 
"computing" a set C, we mean generating increasingly precise "images" of C. At least for 
now, we restrict our attention to bounded subsets of M". 

Consider the two-dimensional case, which is closely related to computer graphics. Intu- 
itively, in this case, the set C is computable if we can draw arbitrarily good "zoom-ins" into 
it. One can view a 2""-precise image of C on a screen as a collection of radius-2~" pixels 
such that the following two conditions are fulfilled: 

1. If a pixel contains a point from C, then it is colored black. This ensured that the entire 
set appears on the screen. 

2. If a pixel is far (say 2~"-far) from C, then it is colored white. This ensures that the 
picture is a faithful image of C. 

We can take the pixels to be balls of radius 2~" with a dyadic center d G D". Formally, 
we say that C is computable, if there is a machine M((i, n) computing a function from the 
family 

r 1 if 5(c/,2-") nc ^ 

f{d,n)=\ if 5(rf,2-2-")nC = (1) 

or 1 otherwise 

On figure^we see some sample values of the function /. It should be noted that the definition 
remains the same if we take square pixels instead of the round ones. It is also unchanged if 
we replace the ratio between the inner and the outer radius to some a > 1 instead of 2. 




Figure 1: Sample values of /. The radius of the inner circle is 2 

One can also define the time complexity Tc{n) as the worst- case running time of a machine 
M{d,n) computing a function from the family ([Q). Low time complexity means that it is 
easy to zoom into the set C. We say that C is poly-time computable if Tciji) < p{n) for 
some polynomial p. 
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So far, one might have been left with the impression that the computabihty definition 
above is not very robust. As seen below, on the contrary, it is equivalent to several other 
reasonable definitions. 

First, we establish a connection between the computabihty of functions and sets. The 
first attempt is to mimic the discrete case, and to say that a set is computable if and only if 
its characteristic function xc is computable. Unfortunately, this is not true: a computable 
function must be continuous, hence xc can only be computable in the trivial cases of C = 
and C = M". 

The next most natural candidate is the distance function 

dc{x) = inf \x — y\. 

y&C 

Unlike the characteristic function, the distance function is continuous and even satisfies the 
Lipschitz condition \dc{x) — dc{y)\ < \x — y\. In fact, we have the following (see [Brv04 j for 
a proof): 

Theorem 2 1. A bounded set C is computable if and only if the distance function dc{x) 
is computable. 

2. If the distance function is poly-time computable, then the set C is poly-time computable. 

3. If n = 1, the converse to part 2 holds: if the set C is poly-time computable, then dc{x) 
is poly-time computable. 

4- If ^ 2, the converse to part 2 holds if and only if P = NP (which is extremely 
unlikely). 

Another possible view on the computabihty of sets it through global computabihty. That 
is, instead of trying to decide one pixel, we are trying to generate an approximation of 
the entire set. The word "approximation" here is in the most natural semi- metric on the 
bounded subsets of R"^, the Hausdorff metric. 

The Hausdorff distance between two bounded sets S and T is the smallest quantity d by 
which we need to "blow-up" 5* to cover T and vice versa. Formally, for bounded 5* and T, 
the Hausdorff distance is 

H{S,T) = inf {d : S C B{T,d) and T C B{S,d)}, 

where B{S,d) = {x : \x — s\ < d for some s G S}. On figure |21 is an illustration of the 
Hausdorff metric. 

On the right-hand image the Hausdorff distance is large, because we would need to blow- 
up A by a large quantity to cover B. As a result, the set _B is a good picture of A on the 
left, but not on the right. 

The theorem below states that C is computable if and only if one can approximate it in 
the Hausdorff metric with finite unions of dyadic balls (see |Brv04j for a proof): 
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Figure 2: An illustration of the Hausdorff metric 



Theorem 3 A bounded set C gMJ^ is computable if and only if there is a machine M{n) that 
on input n computes a finite sequence of centers Ci, 02, ■ ■ ■ ,Ck G lU" and radii ri, r2, . . . , G D 
such that 



In other words, U^^iB{ci,ri) is a good Hausdorff approximation of C. 

It is not hard to see that when time complexity is concerned, we cannot compare the two 
definitions. For example C = [0, 1] x {0} C is obviously poly-time computable, but one 
needs exponentially (0(2'^)) many balls to approximate it within 2~^. 

In order to deal with the computability of bounded sets it suffices to discuss compact (= 
closed and bounded) sets: 

Lemma 4 A bounded set C C M" is computable if and only if its closure, C , is computable. 

We restrict our attention to bounded sets in order for theorem El to make sense. We can 
have the same discussion for unbounded sets by considering the stereographic projection of 
to the n-dimensional sphere. 

2.5 Examples of Computable Sets 

The first natural examples to consider are the simple geometric objects such as a point, a 
line segment, a circle, etc. For these basic examples set computability connects naturally to 
number computability: 
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Claim 5 1. A singleton C = {c} is computable iff all the coordinates of c are computable 
numbers. 

2. A line segment connecting the points x and y is computable iff all the coordinates of x 
and y are computable. 

3. A ball B{c,r) is computable iff r and all the coordinates of c are computable. 

This list can be continued with any other standard shapes. In general, such a shape is 
computable if and only if its parameters are. 

There is another connection between computable functions and sets that gives rise to a 
large family of computable sets. Recall that for a function / : D C M'^ — *• its graph is 
defined by 

We have the following direct connection between the computability of the function / and of 
V f (as a set). 

Theorem 6 Suppose D is a computable closed and bounded set, and f is continuous. Then 
f is computable if and only ifTj is computable. 

A version of this theorem has been first proven in [ZE96] . see also |Brtn8aj and |Brvn4j . 
In section 0] we will give a simple proof of it. 

In particular, the graphs of all the common continuous functions on closed intervals are 
computable. E.g. the graph of f{x) = on the [0, 1] interval. 

Another interesting family of computable sets are the self-similar fractal images. The 
most famous set in this family is probably the Cantor set C. To define the Cantor set let 
Cq = [0, 1] be the [0, 1] interval. Let Ci be the set obtained from Cq by removing its middle 
third: Ci = Co\(|, |) = [0, |] U [|, 1]. We then remove the middle thirds from each of the 
two intervals of Ci to obtain C2 = [0, |] U [|, |] U [|, |] U [|, 1]. Continue this process to 
obtain a chain of closed sets Co D Ci D C2 D . . .. Define C = fl^QCj. See figure El for a 
graphical illustration of the construction. 

The Cantor set has a fractal structure because each of its halves is similar to the entire 
set C with a factor of |. C has an irrational Hausdorff dimension of log3 2, which is smaller 
than 1 but bigger than 0. 

We establish that C is computable by showing it is approximable in the Hausdorff met- 
ric. Ci is a 3~'-precise approximation of C in the Hausdorff metric, and it is very easy to 
approximate Ci in the Hausdorff metric, since Ci is just a union of 2* simple intervals. Thus 
C is easily computable. 

Another famous computable fractal is the Koch snowflake K. The Koch snowfiake is 
obtained from an equilateral triangle by continuously replacing each side of length / by four 
sides of length |, as seen on figure El is a set of Hausdorff dimension logg 4, which is less 
than 2 but more than 1. K is the union of three self-similar sets (corresponding to the sides 
of the original equilateral triangle). 
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Figure 3: The construction of the Cantor set C. 
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Figure 4: The construction of Koch snowflake K. 



As in the case of the Cantor set, the z-th stage of the construction is a 2~®'-*^-approximation 
of K, and it is easy to compute the i-th stage which is a simple union of hne segments. This 
shows that K is Hausdorff approximable, and hence computable. 

An intriguing family of computable quasifractals, the hyperbolic Julia sets is discussed 
in |RW03j . |Ret04j and |Brv04j . Other discussions on the computability of Julia sets, and in 
particular existence of uncomputable Julia sets can be found in |BY04j . |BBY04j . 
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3 The BSS Model 
3.1 The Model 

The BSS model is quite different from the bit model of computability. Like the bit model, 
it also extends the standard Turing machines to deal with the continuous reality. In the bit 
model the extension is through application of the standard machine to continuous problems 
using oracles and naming systems for continuous objects (cf. |Weinnj ). The BSS approach 
extends the model of computation itself. We present here an informal description of the 
model, which is equivalent to a formal one, but is simpler to comprehend for a reader who 
is new to the subject. 

The BSS model in general is formulated for computation over an arbitrary ring or field 
R (for our purposes one can take i? = R or i? = C). The machines in this model are allowed 
to store an entire element of R in one memory cell. The operations the machine is allowed 
to perform on numbers are (1) the ring operations (+, — , ■, and -h if i? is a field); and (2) 
branching conditioned on exact comparisons between numbers (= and <, <, if i? is ordered). 
Initially, the program is allowed to have only some finite number of constants from R. A 
machine computes a function / : i?^ — > i?^ on a domain D C R'^, if on an input x G -D, 
it outputs f{x) E R^ . A machine decides a set C C R^ if it computes the characteristic 
function xc{.x) = x E C, and Xc{x) = otherwise. 

In the BSS model functions and sets of increasing dimension n are sometimes considered. 
The underlying structures then can be viewed over R°° = W^^qR^. The running time is 
allowed to depend on n. In our discussion we restrict attention to the fixed-dimensional 
case, but it can be easily extended to R°°. 

In the case i? = Z2 (the finite field with 2 elements) this machine is equivalent to the 
standard Turing machine. The same statement holds for any finite R. In this sense, the BSS 
model extends the standard computability model. It can be shown that in the case = R, 
the model stays the same if we only allow the machine to have finitely many registers (see 
jB^^.^89j and |Mich89j ). 

In general, all the intermediate computation results of a BSS machine M, as well as the 
output, are rational functions r{x,c) of the inputs x and the constants c of M. From now 
on we restrict our attention to i? = R. We will need the following definition (see 
for more details): 

Definition 7 A semi-algebraic formula 0( finite combination of polynomial 

equalities and inequalities overM. linked by the logical connectives A ("and"), V ("or"), and 
^ ("not"). 

Obviously, any semi-algebraic set is computable by a BSS machine (in constant time). 
As a partial converse, we have the following theorem (cf. Theorem 1 in [ BCSS98j ). 

Theorem 8 If a set C C R" is decided by a BSS machine M , then C is a countable disjoint 
union of semi-algebraic sets. 

A similar statement can be made about BSS computable functions. 
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3.2 Examples of BSS Computable and Uncomputable Sets 



The richest family of examples of BSS computable sets are the semi-algebraic sets. In 
particular, any singleton C — {c}, any line segment and any ball in is BSS computable. 
For example, a ball B{x, r) is given by the simple algebraic condition 

B{x,r) = {ye | (t/i - x^f + (1/2 - + ... + {yn- Xnf < r^}. 

The BSS computable functions include the rational and piecewise rational functions with 
finitely many computable pieces on BSS computable domains. In particular, the step function 
Sq{x) which was not computable in the bit model is easily computable in the BSS model: on 
an input x check whether x > 0, if yes output 1, otherwise output 0. 

In the BSS model, any singleton {c} is computable. In particular, if we take c to be some 
uncomputable number, the singleton {c} C M is BSS computable, but not bit-computable. 
Similarly, the constant function f{x) = c is BSS computable for any c e M, regardless of 
whether there is a feasible way to approximate the values of /. 

We will now present a bit more subtle example of a BSS computable set which is not 
bit-computable. It will be useful later in the discussion. First, it is well known that there 
is a computable (binary) predicate R{x,y) such that the predicate H{x) = 3y R{x,y) is 
uncomputable and if H{x) holds, then the y satisfying R[x,y) is unique. One can think 
of X as the encoding of a Turing machine, and R{x, y) = "the machine encoded by x halts 
after exactly y steps" . Then R{x, y) is computable by a simple simulation, while H{x) is the 
halting problem, well known to be undecidable. 
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Figure 5: A BSS computable set Co one cannot draw. 
We construct the following closed set Co C [0, 1] x [0, 1]. Denote li — [r^^, \] for i — 
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1, 2, . . .. Then [0, 1] = U^Ji U {0}. Define 

Co = ({0}x [0,1])U U hxly. 

R{x,y)=l 

It is not hard to see that Co is closed (there are no accumulation points on the (0, 1] x {0} 
because for each value of x there is at most one value of y such that R{x, y) = 1). See figure 
El for a schematic construction of Co- 

First we observe that Cq is BSS computable using the following program. 

On an input {x,y) G [0, 1] x [0, 1]: 

1. Check whether a; = 0; if yes - output 1. 

2. Otherwise, check whether y = 0; if yes - output 0. 

3. Otherwise, find the rectangle(s) Jj x Ij to which {x,y) belong. This can be done using 
integer numbers and exact comparisons. There are at most 4 such rectangles, if for 
one of them R{i,j) = 1 - output 1, otherwise - output 0. 

Observe that one cannot draw a good image of the set C. In fact, in order to decide 
whether to put a pixel in a small neighborhood of the point (^-^j^, 0^ (this is the middle of 
the interval Jj on the x-axis), one needs essentially to compute H{i), which is impossible. 

On the other hand, there are many functions (and sets) that are uncomputable in the 
BSS model, but are computable in the bit model. For example, the exponential function 
/(x) = on the interval [0, 1] and its graph are uncomputable (see |Brt03bj ). In general, one 
can expect a bit-computable function without some "algebraic structure" to be uncomputable 
in the BSS model. This is caused by the advantage given to the algebraic operations (+, — , 
• and -i-) in the BSS model. 

It is known that BSS computable sets must have an integer Hausdorff dimension (this 
follows from theorem |SI). As a result, the Cantor set (H.d. log3 2) and the Koch snowfiake 
(H.d. log3 4) are not computable in the BSS model. 

3.3 Possible Modifications to the BSS Model 

In this section we discuss possible modifications to the BSS models which address some issues 
raised in the previous section. 

Modifications to the model have been discussed in |B V98j . jBrt98j leading to a notion 
of "feasible real RAM" which is essentially equivalent to the bit-computability. The main 
idea there was, that exact comparisons are not possible on real-life devices, and should not 
be permitted in the model. This idea, by its nature, bars the step function from being 
computable. This is not good when one is concerned with classifying functions into "easy" 
and "hard" , or when discussing computability of multi- valued functions. 
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3.3.1 Uncomputable Constants 

The first concern to address is the use of uncomputable numbers. It is unreasonable to say 
that a function f{x) = a, where a is a constant encoding the halting problem is computable. 
The simple solution is to restrict the BSS machines to use only computable constants. 

Recall that the computable numbers, which we denote by C, are the numbers that can 
be approximated arbitrarily well on a computer. It is not hard to see that C is a countable 
real closed field, and that C + zC C C is algebraically closed. Thus, it makes sense to discuss 
BSS machines over C rather than on the entire M. To emphasize the field C we are working 
with, we will denote this model by BS Sa- 
lt now follows easily that simple geometric objects, such as singletons, line segments and 
balls are -BS^c-computable if and only if they are bit-computable. 

3.3.2 Computation Errors 

The next possible modification addresses the problem of the uncomputability of functions 
such as e^. The BSS model is in part based on the fact that real-life computers usually 
use the four arithmetic operations as a base to performing real computations. However, 

= ^ '^^^ viewed as an infinite-degree polynomial, and is approximated arbitrarily 

well with the finite degree polynomials Pn{x) = Yl^=o^- fact, real-life programs never 
compute e^', but only Pn{x) with some suitably chosen n. 

We further modify BSSc by allowing the machines to err within a given precision e. 
We denote this model by BSS^. In this model, a BSS machine M{x,e) is said to compute 
/(x), if on an input {x,e), e > 0, it outputs f{x) with an error of at most e, and using only 
computable constants. Note that the simple step function so{x) (as well as any other BSSc 
computable functions) is computable in BSSq. 

One can also define naturally the BSS^ computability of sets. A bounded set C is BSS^- 
computable, if there is a BSS machine M{x, e) which uses only computable constants, and 
on input (x, e) it outputs 1 if a; G C and if d{x, C) > e. With this definition, the graph 
of becomes computable. The simple fractals, such as the Cantor set and Koch snowfiake 
mentioned in section also become computable under the modification. 

It should be noted that if we drop the requirement of using only computable numbers, all 
the bounded sets are easily seen to be computable. It is not hard to encode all the "pictures" 
of any bounded set C into one (possibly uncomputable) number c. In other words, all the 
sets are computable in BSS'^. 

3.3.3 Unbounded Computation Branches 

The third modification addresses the problem highlighted by the example Cq on figure|31 The 
problem is the excessive power the BSS model gets from the possibility of having arbitrarily 
long computation paths (as it happened in the example above). 

In the case of "simple" computations, such as e^, or its graph function, we can easily 
estimate the number of steps the machine would have to perform as a function of e. We 
include this condition as an additional restriction on the BSSr machines. 
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We say that a function (or a set) is BSS^"'^ computable, if it is BSS^, computable by 
a machine M, and the running time of M can be bounded by r ( [^J ) for some integer 
computable function r : N — N. We say that t(2"') is the time complexity of the set. 

Under this restriction the set Cq on figure El is not computable, while the function e^' on 
[0, 1], the graph of and the step function so{x) are computable. 

As before, if we remove the restriction of using only computable constants, all the sets 
become computable. 

We summarize the modifications to the BSS model in the following diagram: 

BSSc C BSS'c D BSS'/ 

(Co comp., not comp.) (Cq comp., comp.) (Cq not comp., comp.) 

n n n 

BSS C BSS' D BSS''^ 

{Cq comp., not comp.) (everything is comp.) (everything is comp.) 

We will now show that SS'S'^'^'-computability is equivalent to bit computability for 
bounded sets. Note that they are still different for functions, since the step function So{x) is 
i?5'5'^'''-computable, but not bit-computable. In section El we will connect BSS/ function 
computability to bit computability. 

£ b 

3.4 Computability of sets in BSS^ 
We show the following: 

Theorem 9 Let C dM.^ he a hounded set. Then C is B S S^'' -computahle if and only if it 
is hit-computahle. 

Proof: C is bit-computable ^ C is i^S'S'^'^'-computable 

This is the easier direction. Given a machine M{d, n) for bit-computing C, we show 
how to -BS'S'^'^'-compute it. On an input e, find n such that 2~" < |. Also, by a simple 
binary search, on an input x to the BSS machine, we can find a dyadic c? G ©'^ such that 
|x — (i| < 2~^. We can now use the BSS machine to simulate M{d,n). We claim that the 
output is a valid answer. 

li X e C, then x G 5(c?,2-") n C ^ 0, and M{d,n) outputs 1. If d{x,C) > e, then 
d{d, C) > d{x, C) - \d - x\ > e - 2"" > y > 2 ■ 2~", and M{d, n) outputs in this case. 

Finally, we have to bound the running time as a function of e. This can be done by 
simulating all the possible runs of M{d,n) with some n such that 2^" < |, and all d with 
denominator of 2~("+^+^). 

C is -BS'S'^'^-computable ^ C is bit-computable 

This is a more involved direction. The reduction we will give is not uniform in C. It 
cannot be uniform due to the fact that BSS/ uses arbitrary computable constants. Even 
the simple decision questions, such as equality, are not decidable for arbitrary computable 
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reals presented by Turing machines computing them. Denote the BSS^''' machine computing 
C by M{x,e). 

The nonuniform information needed. Suppose that the BSS machine M uses / con- 
stants ai,a2, ■ ■ ■ ,ai eM. We would need the following algebraic information about ai, . . . ,ai, 
in addition to the Turing machines approximating them: 

1. The algebraic degree of a, over Q(ai, . . . , aj_i), 

2. if this algebraic degree is finite (i.e. Oj is algebraic over Q(ai, . . . , aj_i)), the minimal 
polynomial Pi{x) G Q(ai, . . . , ai_i)[a;] with leading coefficient 1, that has as its root. 
Pi is presented symbolically, with nonleading coefficients given as rational functions 
with non-zero denominators. 

Lemma 10 Provided the nonuniform information as above, for any symbolic polynomial 
p{xi, X2, ■ ■ . , xi) G Q[xi, X2, • • • , x/] we can check whether p{ai, a2, ■ ■ ■ , ai) = 0. 

Proof: We prove the lemma by induction on /. 

Basis: For / = 0, p is just a rational number, and it is trivial to check whether p = 0. 
Step: Assume it is true for / = i — 1, and prove it for I = i. Write 

p(ai, . . . ,ai_i,ai) = pt{ai, . . . ,ai_i)a- + . . . + pi(ai, . . .,ai-i)ai +po(cn, • • • (2) 

There are two cases: 

(z) If Qi is algebraic of degree d over Q(ai, . . . , ai_i), then we can use the minimal poly- 
nomial for ttj to symbolically rewrite (j21) as a degree (at most) d — 1 polynomial in a,: 

p{ai, . . . , ai_i, at) = qd-i{ai, . . . , ai_i)af"^ + . . . + gi(ai, . . . , aj-i)ai + go(ai, • • • , Oi-i)- 

Here the g^'s are rational functions which have non-zero denominators at (ai, 02, ... , flj-i). 
By the minimality of d, p{ai, . . . , ctj-i, Oj) = if and only if all the qj{ai, . . . , aj_i) = 0, 
j = 0,1,. ..,2 — 1, which we can check by the induction hypothesis (since we know that the 
denominators are non-zero, all we have to do is check the numerators at (ai, 02, ... , cti-i))- 
(a) If tti is transcendental over Q(ai, . . . , aj_i), it follows from Q that p(ai, . . . , Oj-i, Oj) = 
if and only if all the Pj{ai, . . . , aj_i) = 0, j = 0, 1, . . . , t, which we can check by the induc- 
tion hypothesis. I 

Lemma 11 Provided the nonuniform information as above, for any symbolic polynomial 
p{xi, . . . ,xi) G Q[xi, . . . ,xi] we can check whether p{ai, 02, ... , a/) > 0. 

Proof: By lemma fTUl we can first check whether p{ai, 02, ... , a;) = 0. If yes, we output 
'no'. Otherwise, using increasingly good approximations, we will eventually be able to tell 
whether p{ai, 02, • • • , o/) > or p{ai, 02, ... , ai) < 0. I 

We now return to the proof of theorem IHl Given a dyadic d G D''' and n G N, we would 
like to compute f{d,n) as in (fT|): 

r 1 if 5(d,2-")nC^0 

f(d,n) = l if 5(d,2-2-") = 

I or 1 otherwise 
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For the rest of the proof set e = 2~". We know there is a computable bound on the running 
time of M{x,e) in terms of e. We compute this bound B = B{e). This means that M(«,£) 
can have at most 2^ different computation paths. Each potential path has an output (either 
or 1), and a set of rational constraints on the input x = (xi, . . . ,Xk) and the constants 
ai, . . . , a/ that ensure that this path is followed. If the constraints are not satisfiable by any 
(xi, . . . , Xk)i it means that the path is never actually followed. The rational constraints can 
be rewritten as polynomial constraints. 

Choose some computation path 7 on which M(», e) outputs 1. We denote the polynomial 
constraints to be satisfied in order to follow 7 by C^{xi, . . . ,Xk,ai, . . . ,ai). We are interested 
whether there is an x G B{d, 2~") fl C. In particular, we would like to know whether there is 
such an x that is accepted by the path 7. This is stated by the following quantified formula 

/^(ai, ...,ai) = 3xi, ...,Xk ((xi - rfi)^ + . . . + (x^ - 4)^ < 2"^") AC^(xi, . . . ,Xfe,ai, . . .,ai). 



Using Tarski's quantifier elimination algorithm (see |Tar51j ). we can convert /^(ai, . . . , a/) 



into a quantifier-free formula (7-^(01, .. . , a/) which has the same truth value. We then can 
use lemmas ITTH and ITD to decide whether g-y{ai, . . . ,ai) is true or false (which is the same as 
deciding /-y(ai, . . . , ai)). A far more efficient procedure can be applied here, using the recent 
developements in algebraic geometry algorithms and the fact that has only some constant 
number of existential quantifiers with no alternations. It is possible to reduce the complexity 
to be exponential in k (and polynomial in the other parameters). See |Ren92j and jBPR03j 
for the algorithms and their analysis. 
As an answer, we output the following 

f{d,n)= y /^(ai,...,a/). 

7 is a 1-valued path of M{»,e) 

As there are at most 2^ such paths 7, the computation will involve computing at most 
2^ times. 

We now show the correctness of f{d, n). First, suppose that there is an x G B{d, e)\^C (re- 
call that e = 2~"). Then M(x, e) must output 1. Let be the computation path correspond- 
ing to X. 7 is an output-1 path, since M(x, e) = 1. Thus, we will have C^(x, ai, . . . , ai) = 1, 
and f-yiai, . . . , ai) = 1. So f{d, n) = 1 in this case. 

Suppose that f{d,n) = 1. It means that there is a computation path of M{x,e) which 
accepts an X G B{d,e). Thus, by definition of SS'S'^'^'-computability, there is a y G C such 
that \x — y\ < e, and d{d,C) < \d — y\ < \d — x\ + \x — y\ < e + e = 2 ■ 2^". Hence 
d{d, C) >2 ■ 2^" implies that f{d, n) = 0, which completes the proof. I 



4 Weak Comput ability of Real Sets 

In this section we introduce another computability and complexity notion for subsets of 
M". This notion was first introduced by Chow and Ko in |CK95j under the name of strong 
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recognizability. It wasn't known at the time that this definition is, in fact, equivalent to the 
standard bit-computabihty definition, as we will show below. 

The idea of the definition is to relax the conditions of the pixel bit computability defini- 
tion. We are given a point x as an oracle to M'^{n), and we must output 1 if x G C and if 
X is 2~"-far from C. Formally, 

Definition 12 We say that a set C is weakly computable if there is an oracle Turing Ma- 
chine M'^in) such that if (p represents a real number x, then the output of M'^{n) is 

( 1 ifxeC 
M'^{n)=l z/S(x,2-") = (3) 

I or 1 otherwise 



Strong Weak 




Figure 6: Strong vs. weak computability. 

Figure El illustrates possible values on input (x, n) for strong (standard) and weak com- 
putability. Note that in the case of weak computability M'^{n) cannot be computing a 
function of x, because there are no nontrivial continuous — 1 valued functions. Hence, for 
some x'es in the 'grey' area the output of M'^(n) depends on the specific oracle for x. 

Our main statement is that, despite the apparent weakness of the weak definition, it is 
equivalent to the standard bit computability. 

Theorem 13 A set is bit- computable if and only if it is weakly computable. 

Before proving theorem ^[ we will give one application of it to obtain a simple proof of 
theorem [HI on the equivalence of function and graph computability: 
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Theorem [Ht Suppose D is a computable closed and bounded set, and / is continuous. 
Then / is computable if and only if is computable. 

Proof: The 'Tj computable ^ / is computable" direction is just a routine graph value 
look-up algorithm and is left to the reader. 

We will prove the opposite direction. Namely, that if / is computable ^ is weakly 
computable, which implies the computability of Fj by theorem 1131 

Given a point (x, y) by an oracle, we run the M'^ computing f onx with precision 2"^"+^). 
We respond to all the oracle queries of M*^ with an answer which is consistent with x as 
well as some point in D. If at any stage of the computation it is impossible, then we've 
discovered that x ^ D, and we output 0. Otherwise, either x & D and we denote x' = x, or 
there is an G -D with |a; — < 2~("+^\ such that the run of the machine on x and x' is 
the same (we don't actually need to know x'). In either case, we obtain a value q such that 
\q - f{x')\ < 2-("+2)_ jf |g _ y| < 2-("+2), output 1; if \q-y\> 2-("+i), output 0. 

To see that the procedure above works, first assume that {x,y) G Fy. Then x & D and 
y = f{x), hence x' = x and \q — y\ = \q — f{x)\ < 2~("+^\ and we will output 1. 

In the opposite direction, if we output 1 then in particular \x' — x| < 2"^"+^'' and \y — 
fix')\ < \y - q\ + \q- fix')\ < 2-("+i) + 2-("+2) = 3 ■ 2-("+2). Hence \{x,y) - {x' J{x'))\ < 
\x' — x\ + \y — f{x')\ < 2~", and d{{x,y),T j) < 2^". Thus we will output whenever 



In terms of time complexity, weak computability is provably weaker than the standard 
computability (assuming P ^ NP). 

Theorem 14 Let C C M" he a hounded set. 

1. If C is poly-time hit computahle, then C is weakly poly-time computahle, 

2. the converse is equivalent to P = NP: "if C is weakly poly-time computahle, then it 
is hit poly-time computahle in general" -vv- 'P = NP", which is most unlikely. Hence 
poly-time hit computahility is strictly stronger than weak poly-time computahility, 

3. a weaker version of the converse holds: if C is weakly poly-time computahle, then it 
is exponential time hit computahle. Moreover, if the machine which weakly computes 
C in poly-time reads at most p{n) hits of the input, then C is computahle in time 

^0(l)20(p(n)+n) ^ 2C(P(")+"). 

See |lBrvn4j for a proof, which we omit here. The last part of the theorem is useful in 
proofs of specific upper bounds, e.g. poly-time computability of hyperbolic Julia sets, see 
|Brvn4j for more details. 

We will now prove theorem IT^ 
Proof: Recall that bit-computability requires computing a function from the family 



d{{x,y),Vf)>2 



— n 




1 if 5(rf,2-") nc ^ 

if 5(rf,2-2-")nC = 



or 1 otherwise 



(4) 
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while for weak computability we need to compute a "function" from the family: 



r 1 ifxeC 
M'^(n) = < if 5(x,2"")nC = (5) 

or 1 otherwise 

C is computable ^ C is weakly computable. 

This is the easy direction in the proof. Suppose we have a Turing Machine computing 
an f{d, n) from the family In order to weakly compute C, we first query d = (f){n + 2), 
and then return /(c/, n + 2). We need to show that (0) is satisfied. 

If a; e C, then \x-d\ < 2-("+2) implies that x e 2-("+2))nC, so E(ci, 2-("+2))nC ^ 0, 
and by (jH) f{d,n + 2) returns 1. 

If B{x, 2-") n C = 0, then \x - d\ < 2-("+2) i^nplies that B{d, 2"" - 2-("+2)) n C = 0. So 
B{d, 2 ■ 2-("+2)) n (7 C B{d, 2-" - 2-("+2)) n C = 0, and by © /(rf, n + 2) returns 0, which 
completes the proof. 

C is weakly computable ^ C is computable. 

We will show the implication for a one-dimensional set C, the proof in > 1 dimensions 
works out in a similar fashion. For convenience purposes assume that C C [|, |], and hence 
in (0)) we only need to consider d G [0, 1]. The proof extends trivially to bigger intervals {C 
is bounded). 

We construct an infinite tree T. In every vertex of T we write a dyadic number. The 
numbers on level / are dyadics of the form m ■ 2~K The root, which is on level 1, is labeled by 
I = 0.1 (all the numbers in this section are in binary notation). Each vertex f on a level / has 
3 children. If the label of f is m ■ 2~' then the labels of its children are m ■ 2~' — 2"^'+^-', m ■ 2~' 
and m-2'^ + 2~'^^+'^\ or in other words (2m - 1) ■ 2-('+i), 2m ■ 2-('+^) and (2m + 1) ■ 2-^^+^\ 
On figure [7| we see the first three levels of the tree (cf fWeiOOj . section 7.2, signed digit 
representation). 

It is easy to see that numbers on every path p in the tree converge to a real number 
Xp G [0, 1]. Conversely, for every x G [0, 1] there is a path p such that Xp = x (e.g. choose p 
to be the prefixes of the binary expansion of x). Moreover, if we denote by p{n) the label of 
the n-th node on p, then \x —p{n)\ < 2^" < 2^'^"-^^\ Hence 0(?t.) = p{n + 1) is a valid oracle 
for Xp. 

We will now describe how to compute f{d,n) as in (jH). On an input {d,n), find two 
nodes vi and V2 on level n such that label(ui) < d < label(t>2) and |label(wi) — d\ < 2~"', 
|label(f2) — d\ < 2~" (if d is an integer multiple of 2~", we can choose Vi = V2)- Denote 
the paths from the root to Vi and f 2 by pi and p2, respectively. Denote the machine weakly 
computing C by M'^{n). We simulate the computation of M'^{n) on the subtrees with roots 
f 1 and V2 as follows. 

Consider the simulation with root vi (the simulation with root V2 is done in the same 
way). For every oracle query 0(m) with m < n, we return the value of pi{m + 1) (which 
is a valid output for the oracle). Otherwise we consider all the possible descendants of Vi 
on level m + 1, and create a separate computation for each of them (thus creating s'^-^+i 
computations). Consider one of the copies and denote the path leading to the selected vertex 
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Figure 7: The first three levels of the tree T 

on level m + 1 by p (p extends pi). If we are now asked about (j){r) for some r < m + 1, 
we return the value of p{r + 1), and otherwise we again consider all possible descendants of 
p(m+ 1) on level r + 1, and split the computation into S*""™ computations. We continue this 
process until all computations terminate. 

If any one of the computations (starting either from Vi or from V2) returns 1, we return 
1. Otherwise return 0. We need to show two things: 

1. The algorithm terminates. 

2. It gives answers that satisfy (jl}. 

We will be using Konig's lemma. 

Konig's Lemma: If a every vertex in a tree has a finite degree, then the tree is infinite if 
and only if it has an infinite branch. 

Suppose that the computation does not terminate. We can view the entire computation 
as a tree where the nodes are the subcomputations described above and a computation Ci is 
the parent of the 3'^ computations it launches. If the entire computation does not terminate, 
then there are two possibilities: either one of the computations C fails to terminate without 
calling to subcomputations, or the tree of all the computations to be performed is an infinite 
tree. 
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In the first case denote the path in T leading to C by p' . Then p' corresponds to a 
dyadic number rf', and also gives an oracle 0' for d' . Note that C is reached and executed 
by simulating M*^ (n). Hence M*^ (n) does not terminate in this case, contradiction. 

In the second case, by Konig's lemma, there must be an infinite branch in the computa- 
tions tree. Denote the branch by Ci, C2, C3, . . .. That is, Ci calls C2, C2 calls C3 etc. Note 
that each Ci works with a path pi of T and pj+i strictly extends Pi for each i, hence the 
infinite sequence of Ci corresponds to an infinite path p in T. The path converges to a real 
number x G [0, 1], and p gives rise to an oracle for x. By the construction, the sequence 
of Ci,C2,C3,... simulates the computation of M'^{n). Hence M'^{n) does not terminate, 
contradiction. This shows that the algorithm terminates. 

We now have to show the correctness of the algorithm. 

Case 1: B{d, 2~^) fl C 7^ 0. In this case, either vi or V2 has a descending path p in T 
which converges to an x G C . Consider the oracle corresponding to this path. One of the 
computation paths of the algorithm will be a simulation of M'^(n), and hence will have to 
output 1. 

Case 2: B{d, 3 ■ 2^") fl C = 0. In this case all points corresponding to descendants of Vi 
and f 2 are at distance at most 2 ■ 2~" from and hence are at least 2~"'-far from C. Hence 
any computation corresponding to simulating M"^ along any of the paths must output 0. 

Note that we are only able to compute a function satisfying a condition weaker than Q. 
Namely, we compute a function / such that 



It is very easy to use / to compute a function that satisfies (jD). Take 

g[d, n) = f{d - 2-("+^), n + 1) V f{d, n + 1) V /(rf + 2-("+i), n + 1) 

If \d-c\< 2-" for some c G C, then either |c/ - c + 2-("+i) | < 2-("+i), or \d - c\ < 2-("+^), 
or \d — c — 2^("+^)| < 2^*^"+^^ So one of the /'s will return 1. On the other hand, if 
B{d, 2 ■ 2-") n C = 0, then B{d - 2-^''+^\ 3 ■ 2-("+^)) n C = 0, 3 ■ 2-("+i)) n C = and 
+ 2-("+i), 3 • 2-("+^)) n C = 0. Hence g returns in this case. This completes the proof. 

■ 

5 Complexity of Real Functions 

In this section we propose a new definition for the computability and complexity of real 
functions, and establish its connections to bit and BSS computability. 

5.1 Computability of Real Functions 

The main idea arises from the equivalence between the function computability and the set 
computability in case of a continuous function, which was established in theorem IHl We have 




1 if s(rf,2-")nc^0 

if 5(rf,3-2-") nC = 



or 1 otherwise 
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a good notion of bit-computability for sets which coincides with .BS'S'^'^'-computabihty. We 
use it to define a computabihty notion for functions: 

Definition 15 We say that a bounded real function f on a bounded domain from a class T 
is graph computable, if its graph Tj is computable as a set. 

Obviously, one has to be careful about the choice of the class JF which extends the class of 
the continuous functions. For example, under this definition, any function / : [0, 1] [0, 1] 
which has a dense graph in [0, 1] x [0, 1] would be computable under this definition. On the 
other hand, one can reasonably interpret this example by saying that all we can know about 
f{x) in this case is that it is some number in [0, 1] (regardless of the precision x is given 
with), and that is exactly what we can read from the graph. 

The definition above extends naturally to multi-valued functions. One possible candidate 
for the class JF is the class of (multi-valued) functions f : D G MJ^ M™' such that / has at 
most finitely many limits at any point of D. In other words, 

|{x} X M'"nr7| < oo (6) 

for aA\ X E D. We denote this class by J-'f. extends the class of the continuous functions. 
Note that both the step function and the square root function on C, viewed as a function 
from to M^, belong to this class. In the case of a function / : ^ R^, all the /'s of 
bounded variation satisfy condition (jH)). The definition will make sense outside the class Tf. 

By theorem definition coincides with the bit computabihty definition in the case of 
continuous functions on closed domains. In the next section we investigate its connection to 
the BSS model, and show that it extends the SS'S'^'^-computable functions. 



5.2 Connections to the BSS model 

In this section we show that all the BSSi^' -computable functions are graph- computable. 
The converse is not true. For example, let Fg be the graph of the exponential function on 
[0,1]. The characteristic function Xr^i^) is easily seen not to be BSS^'^ computable, even 
though it is graph computable. 

Theorem 16 If a function f : M.^ is BSS^^ -computable and bounded on some bounded 

semi- algebraic BSSc-computable domain D (e.g. D = [ai,6i] x [02,^2] x . . . x [0^,6^] with 
computable endpoints), then it is graph computable. 

Proof: We first note that if a function is i?S'5'^'''-computable on D by some machine 
M{z,e), then its graph is i?S'5'^'''-computable. To see this let {x,y) G M'^+^ and £ > be 
given. First of all check whether x E D, and if not - reject. Next, simulate the run of 
M{x, |) to obtain z, \z — f{x)\ < |. If \y — z\ < |, output 1, otherwise output 0. 

li y = f{x), then I2; — 1/| < | and we will output 1. If d{{x, y),T f) > e, then in particular 
\y — f{x)\ > e, so ly — 2;| > \y — f{x) \ — \f{x) — z\ > e — ^ = ^, and the output is 0. It is 
also not hard to see that the running time is bounded by a computable function of e. 

We use theorem IHl to conclude that the graph F/ is bit-computable, so / is graph- 
computable. I 
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5.3 Complexity of Real Functions 

By now we have established that graph computabihty is a useful notion extending both bit 
computability and a natural modified version of the BSS computabihty. Our next goal is 
to give a reasonable definition for graph complexity of real functions. Should we take the 
complexity of computing G/? Or the complexity of weakly computing it? It turns out that 
neither one extends the bit-complexity of / in the continuous case. In fact, we have the 
following theorem. 

Theorem 17 Let f : [0, l]'^ [O5 1]^ be a continuous function, with a polynomial modulus 
of continuity ij,{n) (\f{x) — f{y)\ < whenever \x — y\ < 2~'^^'^^). Denote the following 
properties: 

(a) The graph Tf is poly-time computable as a set. 

(6) f is a poly-time computable function. 

(c) The graph Tj is weakly poly-time computable as a set. 

Then we have the following: 

1. (a) ^ (c) and (b) =^ (c), 

2. in general, (b) =r- (a) implies P = NP, and 

3. (a) ^ (6) ( and also (c) ^ (6) ), implies that integer factoring and other one-way func- 
tions can be done in poly time (this is weaker than P = NP, yet extremely unlikely). 

Proof: The first claim is quite routine and is left to the reader. 

We sketch the proof of the second claim. The goal is to construct a poly-time computable 
function on the [0, 1] interval for which it is NP-hard to strongly compute its graph. The idea 
is to subdivide [0, 1] into the intervals Jj = [j^, j]- Some of them correspond to boolean 
formulas 0^. For such an i we subdivide the middle fifth of /j into 2'^ subintervals, where k is 
the number of variables in 0j. We associate a spike of / to each thruth assignment satisfying 
(f)i. Then it is easy to compute / in poly-time using only one substitution to (pi, but a strong 
query to would allow us to check for the satisfiability of 0j fast. 

We will now prove the third claim, namely that (a) =^ (6) would allow fast factoring 
of integers. The example we construct is for A; = 1, i — 2. It is unclear whether such an 
example can be constructed for k — £ — 1. 

The idea of the proof is as follows. We construct a function / which has the values of 
some one-way function F encoded into it, so that computing / is as hard as computing F. 
We make the graph of / "easy" to draw by the following construction. Let C be a (small) 
interval on which the value F{n) is encoded. We make / fairly dense on intervals CL and 
CR lying on both sides of C. So that its graph in coarse precision looks like two "walls" on 
CL and CR, and there is no need to worry about the value on C to draw the graph. In fince 
precision, on the other hand, the pixels we can query are so small that to answer them on 
C we would only need to know whether F{n) = A; for a few different fc's, which is easy to 
check. 
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We now give a more formal construction. First consider the following function F : N ^ 
{0, 1}*: on an input n, F outputs a canonic encoding of the prime factorization of n. E.g. 
F(6) =< 2,3 >, F(24) =< 2,2,2,3 >, all encoded properly as binary strings. Then F 
satisfies the following properties: 

• F is well-defined on all N, 

• F is one-to-one, 

• the length of F{n) is linear in the length of n: \F{n)\ < c ■ \n\, and 

• given a string s G {0, 1}* and a number n it is easy (poly-time) to verify whether 
F{x) = s. 



1 1 

i+i ' j 



We subdivide each Ij into 5 



Subdivide the [0, 1] interval into the intervals Ij 

intervals: left - Lj, central- left - CLj, central - Cj, central-right - CRj and right - Rj. Let 
k = [logj]. The dividing points of the intervals are dyadic, and the lengths are as follows: 

Lnth{L) ^ Lnth{R), Lnth{CL) = Lnth{C) = Lnth{CR) = 2-('='=+^+3) « Lnth{L). 

Denote 

CL = [xii,Xir], C = [Xir,Xrl] CR = [Xrl, X^r]- 

/ on CL is defined to have the following properties: 

(a) the range of / is in [-2-^ 2^'=] x [-2-^2-^], 

(b) f{xu) = f{xir) = {0,0), 

(c) for any {x,y) such that \y\ < a < there is a 2; G CL such that \z — x^] < 
a ■ 2-('=^+3) and \ f{z) -{x,y)\<a- 2-(^fc+'=+3)^ 

{d) f is poly-time graph computable on this interval. 

This can be achieved, for example, by a function presented on figure |H1 (left). 
Next, we define / on C as follows: 



(2'^*^ ■ {z — xir), F{n) ■ {z — xir)), if z is in the left half of C 
{2^^ ■ {xri — z), F{n) ■ {xri — z)), if z is in the right half of C 



See figure IHl (right). Note that f{xir) = f{xri) = (0,0). We define / on CR symmetrically 
to its definition on CL above. / is defined to be (0, 0) on the intervals L and R. It is easy 
to see now that / extends continuously to the closed interval [0, 1] by setting /(O) = (0, 0). 

If / were poly-time computable, then querying its value on C for the interval /„ with 
precision 2"*^^'°^"'^ would allow us to compute F{n). Our goal now is to show that the graph 
Fj- is poly-time computable as a set. The idea is as follows: if we query a fairly large pixel 
of Fj, then the picture will be dominated by (the computationally easy) graph of / on CL 
and CR, and we will not have to worry about the value of / on C. If we query a fairly small 
pixel in C, then to determine the answer we will have to verify at most one potential value of 
F, which can be done in poly-time. On an input {d,m), d = {z, {x,y)) G D^, we will decide 
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Figure 8: The construction of f 



whether B{d, 2"™) n T/ 7^ or B{d, 8 • 2-^") n = 0, with the other cases being the "grey 
area" . It is not hard to see that this is equivalent to the standard definition of computing a 
set. 

First of all we find the interval /„ which is relevant to the particular d. If there is more 
than one relevant /„ then the problem easily decidable in poly-time. Next, it is easy to 
decide whether d is close to the graph of / on the intervals L, CL, CR and R. The hard part 
is to deal with the interval C. Since the picture is symmetric with respect to the middle of 
C, we assume without loss of generality that z is to the left from the middle of C. 

U z < xir — 2 ■ 2~"^, then the pixel in question is far from C, and the value of / on C is 
irrelevant. Otherwise, there are two cases: 

Case 1: - 2 ■ 2'"^ < ^ < x;^ + 3 ■ 2""^. 

In this case, in for the value of / on C to become relevant we must have \x\, \y\ < 
2'^^-4-2~'^-|-2~'" and \x\, \y\ < 2"^^. By property (c) of / on CL there is a point p = {z', f{z')) 
such that \z' - xir\ < 2"™ and \ f{z') - {x,y)\ < 2"™, thus |(/, f{z')) - rf| < 5 • 2"™, and it 
is safe to output 1 in this case regardless of the value on C . 

Case 2: z> xir + ?>- 2-"^. 

In this case we only need to care about potential points f{z') on the graph with z' > 
z — 2 ■ 2"™ > xir + 2""*. It is not hard to see that for these values the distance between any 
two potential points on the graph of / (for two different possible values of F{n)) is at least 
2~™/-\/2, hence {x,y) can be only close to at most 4 potential graphs of /, and to confirm 
those we need to query F{n) = tl for 4 t's, which is done in poly-time. This shows that Fj 
is poly-time computable. I 

What we need is a complexity definition that extends part (6) in theorem [T7I to all the 
graph-computable functions. 
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Definition 18 We say that a bounded multi-valued function / : D C M*^ — for some 
bounded computable D is graph-computable in time T{n,m) if there is an oracle Turing 
machine M't'{m,d,n) which, given an oracle for x & D, computes a function from the 
family 

{1, if \d — y\ < 2"" for some y G /(x) 

0, zf\d-y\> 2-2-'' for ally e f iB{x,2-"')) (7) 

or 1 otherwise 

and the computation time is bounded by T{m,n). One can think about M'f'{m, •ju) as 
(strongly) computing a set Am C with precision 2~^ such that 

/(x)cA„C/(i?(x,2— )) (8) 

The set can be thought of as a "vertical cross-section" of the graph of /. In the case 
of a continuous single- valued /, Am {/(^)} as m — > oo. More generally, one can easily 
see that ^ is equivalent to the condition 

oo 

fl / 2-^)) cAmCf {B{x, 2-^) . (9) 

i=0 

For definition to make sense, we need it to be equivalent to definition as far as 
computability is concerned. 

Theorem 19 A function f : D G ^ for some closed and bounded computable D is 
computable as per definitions^ if and only if its graph is computable (that is, it is computable 
as per definitions^. 

Proof: / is computable as per definition ^1 ^ is computable. 

We show that in this case F/ is weakly computable. Given a point {x,y) by an oracle, 
we run the M'^{n + 2, •, n + 2) computing f on x with precision 2~("+^). We respond to all 
the oracle queries of M'^ with an answer which is consistent with x as well as some point in 
D. If at any stage of the computation it is impossible, then we've discovered that x ^ 
and we output 0. Otherwise, either x E D and we denote x' = x, or there is an x' E D 
with \x — x'\ < 2^("+^), such that the run of the machine on x and x' is the same (we don't 
actually need to know x'). In either case, we obtain a set An+2 satisfying condition (jH)) (with 
n + 2 instead of m and n). If d{An+2, y) < 2"("+^\ output 1; if d{An+2,y) > 2"('"'+^\ output 
0. 

If a; G -D and y G f{x), then y G An+2 by (jH)), so d{An+2, y) = 0, and we output 1. 

If we output 1, then d{An+2, y) < 2"*^"+^^ thus by (jHl) there is a 2; G D such that \x' — z\ < 
2-("+2), and \y-f{z)\ < 2-("+i). Thus \{x,y) - {zj{z))\ <\x~z\ + \y-f{z)\ < \x - x'\ + 
\x'~z\ + \y- f{z)\ < 2~("+2) + 2-("+2)+2-("+i) =2"". Hence in the case rf((x, y), F/) > 2"", 
we will output 0. 

Fj is computable ^ / is computable as per definition ITHl 
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Note that by theorem El the reduction in this direction cannot be too efficient. We 
denote by Hi the projection vr^ : R'^ x ^ M^, n£{x, y) = y. 

Given a point x G -D and a triple (m, n), (i G we would like to check whether d is 
2^^-close to an Am as in (jH)). To do this we take a dyadic approximation dx G D*^ of a; such 
that 14 - a;| < 2-(™+2). If d is 5 ■ 2-("+2)_ciose to ne {Vf n 5(4, 2-(™+2)) ^ R'^) we output 1, 
if it is 7 ■ 2-("+2)-far from vr^ {Vf f] B{dx, 2-(""+i)) x R^) we output 0. This can be done by 
"drawing" a portion of the graph using queries for Tf. 

If \d — f{x)\ < 2~" for one of the values of /(x) (in which case we must output 1), then 
/(x) G {Tf n 5(4, 2-(™+2)) x R^), and \d - f{x)\ < 5 ■ 2-("+2), hence we output 1. 

If we output 1, then there is a z G 5(4, 2-('"+^)) such that \f{z) - d\ < 7 ■ 2~("+2)_ 
We have \z - x\ < \z - 4| + |4 - x| < 3 ■ 2-('^+2)_ jj^^g ^f(^^^ _ ^| < 2 ■ 2-^^ for some 
z G 5(x, 2^™"), and we are allowed to output 1 in this case. I 

We will later see that definition ^1 extends the standard complexity definition for the 
case of the continuous single-valued functions. We start off with some examples from the 
discontinuous/multi- valued functions. 



5.4 Examples of Function Complexity 

We begin with the simplest discontinuous function, the step function so{x) = if x < 0, and 
so{x) = 1 if X > 0. It is now easy to see that sq is graph-computable in linear time (under 
definition ITHj). We ffist query for x with precision 2"^""*"^^. If |x| > 2~^"'~^^^ we output either 
{0} or {!}, depending on the sign of x. If |x| < 2~", we output {0, 1}, so in either case (jH)) 
is satisfied. 

In general, if we have a piecewise continuous single- variable function / with finitely many 
pieces go, gi, . . . ,gk and computable discontinuity points oi, a2, . . . , a^, then / is computable, 
and the time complexity T(/) satisfies 

k k 

T{f)<Y,T{gd+Y.^{a,)+0{n). 

i=0 i=l 

Similar relations can be derived for functions over higher dimensions. 

Consider the square root function : C — C. It is two- valued at all points but 0. It 
cannot be made computable on C in the bit-sense because it does not have a continuous 
branch defined on the entire C. It is also uncomputable in the original BSS model, but has 
an efficiently computable branch in BSSq'^. 

The following algorithm graph-computes the square root in poly-time on a bounded 
domain: 

On an input x, (m, d, n): 

1. Check whether |x| < 2^*^^"+^^ In this case /(x) with precision 2~"' looks like the point 
{0}. 

2. Otherwise, we can compute r and such that x = re^'^*'^ with precision 2"^^"'+"^) in 
poly-time. 
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3. compute r' = y/r, the positive real root, 



4. take one of the values of 0' = | (mod 1), 

5. we take = {r'e^'^^'^', -r'e^'^*'^'}. 

It is not hard to see that other simple multi-valued functions, such as x i— > ^ on the 
complex plane are computable in poly-time. 

Next, we consider the characteristic function xa for a bounded A C M". Note that the 
computability of xa is equivalent to the computability of F^^, which is equivalent to the 
computability of A and A'^. 

Consider the complexity of xa- On an input x we need to output one of the three possible 
sets: {0}, {1} or {0, 1}. We must include 1 ii x E A, and must exclude it if x is 2~'^-far from 
A. Similarly, we must include ii x ^ A, and must exclude it if x is 2~'^-far from A'^. We see 
that the complexity of xa in this case is roughly equal to the sum of the weak complexities 
of A and A''. 

5.5 The Definition Extends Standard Function Complexity 

In this section we show that for continuous functions with a reasonably small modulus of 
continuity graph-complexity extends the standard complexity definition. In particular such 
a function is poly-time computable if and only if it is poly-time computable according to 
definition ^| We prove the following theorem: 

Theorem 20 Let f : D ^ M} be a continuous function, where D C M.^ is bounded. Then 
the following holds: 

1. If f is computable in time T{n) according to the standard definition, then it is com- 
putable in time T{n + 2) + 0{n) according to definition\l^ 

2. If f is computable in time S{m, n) according to definition \1^ and the modulus of 
continuity for f is a computable /i = yu(n) (so that \f{x) — f{y)\ < whenever 
\x — y\ < 2^^*^")^, then f is computable in time 0{n ■ S{fj,{n + 2),n + 2)) according to 
the standard definition. 

Proof: (1) For an input {m,d,n), compute an approximation q of f{x) with precision 
2-("+2). If \d-q\ < 5-2-("+2), output 1. If \d-q\> 7-2-("+2) output 0. It is not hard to see 
that the result of such computation satisfies ((Tj). The time required here is T{n + 2)+0{n) 
for the comparisons. 

(2) The set A^(n+2) C f{B{x,fi{n + 2))) is contained in 5(/(a;), 2-("+2))^ ^^^^ 

we can 

approximate f{x) within 2~" by 0{n) queries about the distance of d from this set in different 
scales. This is done in time 0{n ■ S{fi{n + 2), n + 2)). I 
In particular, theorem 1201 implies that a continuous function with a polynomial modulus 
of continuity is poly-time computable if and only if it is poly-time computable according to 
definition UHl 
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5.6 Different Bit-Complexity Notions for Discontinuous Functions 

In this section we compare the three different complexity notions for general (possibly dis- 
continuous and multi- valued) functions. In theorem 1171 we have seen three poly-time com- 
putability notions, in the general case they correspond to the following: 

(a) The graph F/ is a poly-time computable as a set. 

(b) / is a poly-time computable function as per definition 

(c) The graph F/ is weakly poly-time computable 

In the continuous case we have seen that (a) =^ (6) =^ (c), while the converses are 
extremely unlikely. In the general case we have: 

Theorem 21 /n the case of a general f : 

1. (a) =^ (c) and (b) =^ (c), 

2. there is a function that satisfies (6) and (c), but does not satisfy (a) unless P = NP, 

3. there is a function that satisfies (a) and (a), but does not satisfy (b) unless P = NP 
(FACTORING ^ P in the single-valued case). 

Proof: The first two parts are standard constructions, and are left to the reader. We will 
only prove the third part. 

The proof uses the same idea as the construction in theorem [T71 but is much simplier due 
to the fact that we are now allowed multi-valued, discontinuous functions. The multi-valued 
function we construct is / : [0, 1] — > P{[0, 1]). 

For each j eN set f{j) = 

[0,1], also set /(O) = [0,1]. 

Let n be some natural number. If n does not encode a boolean formula, set f{x) = {0} 
for all X G [^^, ^] • Otherwise, n encodes a boolean formula where x is a vector of 

k < log n — 1 variables. 

We create a 1 — 1 correspondence between points of the form | + G [|, l], t = 
0, 1, . . . , 2^^ — 1 and the possible assignments for (p. For any x G [;^, ^] , set 

f{x) = {0} U I - + 7-7—7 : t corresponds to a satisfying assignment for 

A graph of / is schematically presented on figure El 

We show that Fj is (strongly) poly-time computable. Let (^1,^2) G and n be given. 
We need to decide whether (^1,^2) is 2~"-close or 2 ■ 2~"-far from (6/1,^2). If di ^ [0,1], 
or d2 ^ [0, 1] it is very easy to answer the query. Suppose di G [^^, ^] for some m. If 
^ < 2"" we can just output 1, and we do not need to compute /. Otherwise, we check 
whether m corresponds to a boolean formula in variables. If it doesn't, then it is very 
easy to answer the query. Otherwise, since 2~" < ^ ^ 2^^' ^^^^ need to make at most 
two substitutions to in order to correctly answer the query. 

On the other hand, if / were graph poly-time computable as in (6), we would be able to 
decide satisfiability in poly-time, which would imply that P = NP. I 



Subdivide the interval [0, 1] into the intervals / 
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Figure 9: The construction of f 



5.7 Notes on Complexity in the BSS Model 

In theorem 1161 we have seen that from the computability point of view, graph computabihty 
extends BSS^j' -computabihty of functions. In this section we give an example that illustrates 
the gap between the definition from the complexity point of view. 

We give an example of a set A C [0, 1] x [—1, 1] which is BSS computable by a machine 
M without constants or error, such that the running time of M on (x, y), < a; < 1 is 
bounded by a polynomial in — logx. Thus the function 

f{x,y) = X ■ XAix,y) 

is BSSq^ computable, with a time bound polynomial in — \oge (this is the natural "poly- 
time computability" for this model). On the other hand, graph-computing / in poly-time 
would allow us to solve SAT in P, and imply that P = NP. 

The set A is defined as follows: it is a union of Aj's which are small sets in the neigh- 
borhood of (7,0), and i = 2s, where s is a number representing a boolean formula with 
k < log 2 variables. We associate the 2'' assignments to with numbers from to 2^^ — 1. Ai 
is defined as follows: 



A, 



11 V 

{x,y) : - < a; < - — -, y>0, < — 

t t — 1 X — l/i 



< 1, and 



2^-y 



X — l/i 

On figure ITUl is an illustration of A^ corresponding to X2, X3) = (xi A X2) V X3. 
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Figure 10: Ai corresponding to (f){xi, X2, X3) — {xi A X2) V xs 

To compute A in the BSS model one first locates the i such that ^ < x < j^. This takes 
time linear in — logx. Then it takes time polynomial in logi ~ —logx to verify whether 
{x,y) G Ai. Thus A is BSS computable in time polynomial in — logx. 

Now suppose that f{x,y) is poly-time computable in the graph model. Then a run on 
X — j^, y = with 2""* < 2^ will reveal whether or not the formula (j) corresponding to s is 
satisfiable, yielding a poly-time algorithm for SAT. 
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